AN UPPER BOUND ON THE HIGGS BOSON MASS 
FROM YUKAWA UNIFICATION 
AND A COMMENT ON VACUUM STABILITY 

CONSTRAINTS 

Nir Polonsky * 

Department of Physics, University of Pennsylvania, Philadelphia, Pennsylvania, 19104 

^' 

0\ ■ May 1994 

0\ ' 



Abstract 



Only small regions in the mt — tan /3 plane are allowed when considering simultaneously 
(assuming the MSSM) coupling constant unification and (minimal) GUT relations among 
^ ' Yukawa couplings (i.e., ht = hr at the unification point). In particular, if mt ^ 175 GeV 

i we find that only 1 ^ tan/3 ^ 1.5 or tan/3 ^ 40 ± 10 is allowed. The former implies that 

the light Higgs boson is ^ 110 GeV and, in principle, visible to LEPII. The prediction for 
the Higgs boson mass in the tan /3 « 1 scenario is discussed and uncertainties related to (i) 
' vacuum stability constraints, (ii) different methods for calculating the Higgs boson mass, 

T^lj- . (in) two-loop calculations and (iv) GUT corrections are briefly reviewed. It is shown that 

OS ' large left-right mixing between the t-scalars can significantly enhance the Higgs boson mass. 

That and an ambiguity in the size of the two-loop correction lead to our conservative upper 
pi i' bound of 110 GeV. Vacuum stability considerations constrain the t-scalar mixing and slightly 

JL . diminish the upper bound (depending on the value of mt). Improved two- loop calculations 

Qj are also expected to strengthen the bound. 

Realizing the minimal supersymmetric extention of the standard model (MSSM) within a 
grand-unified theory (GUT) and assuming minimal matter content, often imply specific relations 
^ I between the different Yukawa couplings. When considering the simplest example, i.e., "bottom- 
ed ■ tau unification" [Q] hf, = h^-, we find Q (in agreement with others [^]) that the predicted ranges 
for as{Mz) and for tan/? = Vup/i^down are strongly correlated with the weak angle s^(M^) pa 
0.23240 and with mt, and are strongly constrained]^. Such relations provide a consistency test 
of simple GUT models, as well as a means to probe the GUT structure. Furthermore, they can 
be used to distinguish such models from generic string models that typically do not constrain 
the Yukawa couplings in this manner. Thus, exploration of all implications of the bottom-tau 
unification (and similar relations) is well motivated. Below, we study the implications for the 
mass of the light MSSM (CP even) Higgs boson m^o. We will motivate the exploration of the 
tan /3 — > 1 limit and show that m^o in this region is heavier than naive expectations, but most 
probably not too heavy to be seen in LEPII Exploring that region, where the SM minimum 

'Talk presented at SUSY- 94, International Workshop on Supersymmetry and Unification of Fundamental 
Interactions, Ann Arbor, Michigan, May 14 -17, 1994. Pennsylvania Report No. UPR-0595T. 
iPor mf'" ^ 143 GeV. 

^ Us, tan/3 and mt are the strong coupling, the ratio of the two Higgs doublet expectation values and t-quark 
mass, respectively. 
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oc cos^ 2/3 is shallow, requires us to pay special attention to vacuum stability constraints and we 
will briefly comment on that issue. 

The most impressive evidence for TeV scale supersymmetry is the well publicized coupling 
constant unification ||5|. However, it is important to note that [up to matching functions induced 
by GUT-scale physics, which dominate the 0.008 theoretical error bar in (|I|)] unification requires 
as{Mz) ^ 0.12, i.e., 

as{Mz) = 0.125 ± 0.001 ± 0.008 + Ha, + 3.2 x 10-'^GeY-^[{m^°^y - (143 GeV)^]. (1) 

The heavier the t-quark, the more the extracted value for the weak-angle is diminished and the 
predicted value for as is increased. Two- loop contributions from Yukawa interactions Ha, can 
diminish by as much as 0.003 (~ 2.5%) if the Yukawa coupling is of order unity (i.e., the 
function Ha, ^ —0.003). The above observations affect the region in parameter space which is 
consistent with the observational range for the 5-quark mass and with bottom-tau unification. In 
particular, for m\°^^ ^ 200 GeV only two regions are allowed: {%) tan/3 ~ 1 (see Fig. 1) and (ii) 
large tan fj. For example, if mf°'^ ~ 174 GeV then tan /3 « 1.5 or tan /3 « 40ib 10. Hereafter, I 
will discuss region [i) only. An important implication of Yukawa unification, choosing the tan (3 
near unity region, is that the SM-like Higgs boson mass nearly vanishes at tree-level. Thus, a 
careful examination of the loop corrections is required. 

There are important theortical uncertainties which constrain the predictive power of the 
model Q. These are taken into account (assuming no conspiracies), e.g., in (^) and Fig. 1, but 
will not be discussed here. The following points, however, should be stressed: 

• Us dependence: We use the predicted value of as{Mz)[s^ ,mt, ...], e.g., as{Mz) ~ 0.12 
for m^°^'^ « 100 GeV and as{Mz) ~ 0.13 for mf°'^ 180 GeV. Using instead a fixed 
value a^{Mz) for the whole mt range implies some hidden assumptions on GUT-scale 
corrections. The lower the value of ag (and the heavier the t-quark) the larger and less 
likely are these corrections. 

• Finite-loop corrections to mb{MsusY)- The relevance of those corrections was pointed out 
[0] recently. However, while being extremely important for large tan /3 (maybe even to the 
point of removing any predictive power), they are negligible in region (i), as can be seen 
in Fig. 2. The details of the low-energy spectrum are nearly irrelevant in determining the 
size of that region. 

Different aspects of region (i) were discussed by various authors. The possibility that the 
Higgs boson mass is induced only at the loop-level was discussed by Diaz and Haber |p. The con- 
sistency of that scenario with bottom-tau unification recently led to greater attention to that re- 
gion of parameter space ^, |l^. The consistency with Yukawa unification is easily understood: 
From perturbativity one has a mt-dependent lower bound on tan /3 ^ 1. Thus, we are in a region 
of large top- Yukawa coupling ht ~ 1. Slightly decreasing tan/3 (i.e., decreasing sin /3 or increasing 
ht) at Mz would lead to divergences (i.e., ht "hits" its Landau pole) at higher scales [|ll|: ht has 
a quasi-fixed point at tan/3 ~ 1. This balances too large corrections to the hb/hr ratio, can 
explain the heavy t-quark, and makes region (z) relatively insensitive to theoretical uncertainties. 

For tan/3 = 1 there is a custodial symmetry SU{2)l x SU{2)ji — > SU{2)L^ji in the Higgs 
sector p^]. The deviation of tan/3 from unity measures the symmetry breaking, which is induced 
at the loop level due to the Yukawa interactions. Thus, the enhanced symmetry can be respon- 
sible for the smallness of the region. It also puts an upper bound on mf ''^ ^ 185 GeV so that 
tan/3 < 2. 
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As mentioned above, we will not discuss all possible signatures, but will focus on the loop- 
induced Higgs boson mass. (Another interesting signature is a possible light mixed t-scalar 



[|lO|].) In particular, Barger et al. argued |13], using the leading-logarithm approximation, that 
ruf^o ^ 85 GeV for mf°'^ ^ 160 GeV. We will show that this is rather a typical mass range 
and that the upper bound, though still relevant for LEPII, is much higher. An important 
enhancement comes from the large left-right mixing in the t-scalar sector and from the large 
split between the two t-scalar mass eigenstates. That enhancement leads to the upper bound 
nif^o ^ 100 (110) GeV for m\°^^ ^ 160 (175) GeV. The enhancement is sensitive to the type of 
vacuum stability constraints one imposes (see below). 

From the minimization of the Higgs scalar potential one has for the supersymmetric Higgsino 
mass parameter /^^ oc l/[tan^ 1]) and diverges as tan /? — > 1, i.e., in the symmetric limit. (In 
fact, one expects finite loop corrections to be relevant near the limit, and we can assume, without 
loss of generality, tan/3 ^1.1 p^].) The large //-parameter dictates the phenomenology of the 
scenario. In particular, the tree-level CP even mass matrix is nearly degenerate: it has a very 
heavy and a nearly zero mass eigenvalue (this is a trivial consequence of the custodial symmetry) . 
In practice, the light tree-level eigenvalue m^o < Afz| cos 2/?| is small but grows with mf"'^ (i-e., 
with the lower bound on tan/3). At the one-loop level, one has m^o = rrv^o^ + A?o A?o- The 



upper bound on the loop correction can be estimated |l5|| (assuming only t-scalar contributions): 
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Due to the large parameter the left-right mixing O-t between the two t-scalars can be substantial 
and the leading-logarithm and the mixing Ag^ terms in (^) can be equivalent. Thus, mj^o 
is enhanced by an additional factor of ~ \/2. However, since tan/3 increases with mt (from 
perturbativity), |/i| and Ag^ decrease and we obtain an interesting interplay between the overall 
factor of and Ag^. 

Eq. (^) is, of course, only an approximation. The results of a complete calculation for the 
Higgs boson mass are shown in Fig. 3 (for the details of the calculation, see Ref. 4). tan/3 
is constrained as in Fig. 1, we assume the gaussian distribution mf°'*^ = 143 it 18 GeV (from 
precision data^ and all superpartner masses are constrained to be ^ 1 TeV. The monte-carlo 
generated histograms contain information on the upper bounds and on the distributions. The 
upper bound is a function of m^: 

• m^o ^ 100 GeV for mf^ ^ 160 GeV, 

• m^o ^ 110 GeV for mf ^ 175 GeV, 

• m^o ^ 120 GeV for mf''' ^ 185 GeV. 



^Using the recent CDF range[^ for mf"'*^ would only change the relative population in the different histograms 
and would not affect our discussion. 
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The distribution has two peaks, one which is enhanced by the i-scalar mixing term and which 
determines the upper bound, and a peak at a much hghter mass from points with no enhance- 
ment. Thus, the mass is typicahy much hghter than the upper bound, and, if the t-quark is 
hghter than ~ 170 GeV, may stih be relevant for LEPI. (Note that we do not attempt to give 
rigorous bounds.) 

We compare the above bounds to those derived for any tan/3 in Fig. 4. The ~ 30 GeV 
difference is due to the non- vanishing tree-level mass ^ Mz in the general case. (Loop corrections 
are, however, typically smaller in the general case.) One can also compare either case with typical 
upper bounds mf^o ^ 130 — 150 GeV derived from perturbativity considerations [^], and to the 
lower bound given by Sher |17] for the non-supersymmetric case, m/jO ^ 132 GeV. However, 
strong assumptions regarding vacuum stability are made in the latter case. 

Complications in the calculations, which are explored elsewhere Q, suggest even stronger 
upper bounds: 

1. The magnitude of the two- loop terms is known to be small and probably negative 1 18]. 
However, the exact magnitude is ambiguous given the complicated low-energy structure of the 
model. For example. In Fig. 3 we used the effective potential method|19] to extract mf^o. In 
that method A^jO is given to a one-loop leading-logarithm order. It is straightforward to show 
that the renormalization-point dependence of the one- loop leading logarithm [e.g., in (§)] is 
equivalent to two-loop next-to-leading logarithms which are typically positive. Thus, one could 
overestimate when choosing the renormalization point Q = Mz (as in Fig. 3). In Fig. 5a 
Q = 600 GeV, and indeed lower masses are suggested. In particular, if the scalar quarks are of 
the order of TeV or heavier, one would falsely get large Higgs boson masses using that method, 
unless, Q is properly adjusted (so that the loop expansion does not break down). Alternatively, 
we could use different methods, e.g., the renormalization group method [^] leads to typically 
lower one- loop upper bounds (in agreement with Fig. 5a). This is because one is implicitly 
including negative two-loop leading logarithm contributions in that method. In short. Fig. 3 
corresponds to a conservative estimate of the upper bounds. Once the non-trivial task of a 
complete two-loop calculation is carried out, we expect slightly stronger upper bounds. 

2. In some cases large i-scalar mixing corresponds to the SM model populating only a local 
minimum of the full scalar potential. That is, when considering all scalar fields on equal footing, 
the global minimum may be in a direction in which the t-scalar has a vacuum expectation value 
and, thus, the physical vacuum is not stable. Eliminating cases in which a global color and 
charge breaking minimum exists slightly reduces the upper bound and also affects the mixing 
peak population. The effect is reduced with increasing nit and never diminishes the mixing 
peak completely. One could try to identify dangerous directions in the multidimensional scalar- 
field space and characterize those directions by analytic constraints. However, such constraints 
entail assumptions that are not always useful for identifying the deepest minimum (i.e., the 
constraints are too weak). In addition, analytic constraints typically do not distinguish a local 
minimum from a global one and can exclude legitimate points (i.e., are too strong). Thus, 
analytic treatments are often not satisfactory. Furthermore, a point in parameter space which 
is consistent with electroweak symmetry breaking contains at least a local minimum which does 
not respect color and charge symmetries. This is shown in Ref. 4 where we also compare 
numerical and analytic treatments of the problem. Recent progress was also reported by Carena 
and Wagner [21|. 

Lastly, the prediction of m^^o is relatively stable when considering model-dependent evolution 
of the soft parameters between the Planck and GUT scales 1 22 ] . This is shown in Fig. 5 and is 
due to the cancellations between corrections to the relevant soft parameters and to /i^. (However, 
the t-scalar mixing and mass are affected.) Thus, our discussion is only slightly affected by that 
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uncertainty. 

To conclude, we studied the one-loop Higgs boson mass in the region of parameter space 
which is well motivated by Yukawa unification. We treat all theoretical uncertainties in a 
reasonable manner and the upper bounds given, e.g., m^o ^ 110 GeV (for mf"'*^ ^ 175 GeV), 
correspond to the most conservative estimate. The Higgs boson in this scenario is most likely 
in a range that is, in principle, visible to LEPII [unless the t-quark is heavier than ~ 180 GeV, 
in which case there is not much motivation to consider region (i)]. Two-loop calculations and 
vacuum stability constraints slightly strengthen the upper bounds. 
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Figure 1: The allowed region for tan/3 ^ 1. The t-quark pole mass ° ^ is in GeV. The 
lower bound on tan (3 is from ht < 3 at the GUT scale, and the upper bound from bottom- 
tau unification (including theoretical uncertainties). The one s.d. range for m^"^^ , suggested 
by precision data, and tan (3 = 1 are indicated for comparison. Recent CDF analysis suggests 
mf'^ 174 lb 16 GeV. For tan /? < 1 no model with universal boundary conditions is possible. 



Figure 2: The correction to the 6-quark MS mass mb{Mz) ~ 3 GeV (in GeV) from finite 
superpartner-loops. m^°^'^ = 174 GcV and universal boundary conditions for the soft param- 
eters at the unification scale are choosen at random. All points included are consistent with 
electroweak breaking, tan/3 = 20 is indicated for comparison. For tan/3 < 20 corrections are 
less than 10%. 



Figure 3: The Higgs boson mass mj^a (in GeV) distribution in a sample of monte-carlo calcula- 
tions with random universal boundary conditions at the unification scale and using the effective 
potential method with subtraction scale Mz- tan/3 is constrained as in Fig. 1. m^°^^ is in the 
range (a) [155, 165], (6) [165, 175], (c) [175, 185] GeV and has the gaussian distribution 143 ± 18 
GeV. 



Figure 4: The upper bound on the Higgs boson (MS) mass m/jO as a function of ^ (all 
masses are in GeV) . In the constrained case tan /3 is constrained by bottom-tau unification as in 
Fig. 1. The general case is for all tan/3 (consistent with electroweak breaking). Note the local 
minimum in the constrained curve which is due to the interplay between the mixing enhancement 
and the factor. The upper bound curves are derived using monte-carlo methods and are not 
rigorous. Generic upper bounds in the MSSM and in non-minimal models (NMSSM) derived 
from perturbativity considerations are shown for comparison. Also shown a suggestive lower 
bound in the non supersymmetric case (we use ag = 0.125 in the formula of Ref. 17). 



Figure 5: The Higgs boson {rriyfi) vs. Higgsino (/i) mass plane (in GeV) for m^°^ = 174 GeV, 
tan/3 = 1.52, and with random universal boundary conditions at (o) the GUT scale 10^^ GeV, 
(6) at the (reduced) Planck scale 10^^ GeV. In case (6) we assume a minimal SU5 model above 
the GUT scale (for details, see Ref. 22). The effective potential method with a subtraction scale 
Q = 600 GeV is employed, which leads to lower values of m^o than for Q = Mz- 
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Department of Physics, University of Pennsylvania, Philadelphia, Pennsylvania, 19104 

ABSTRACT 

Only small regions in the nit — tan (3 plane are allowed when considering simul- 
taneously (assuming the MSSM) coupling constant unification and (minimal) GUT 
relations among Yukawa couplings (i.e., hf, = h^- at the unification point). In partic- 
ular, if rrit ^ 175 GeV we find that only 1 ^ tan [i ^ 1.5 or tan [i ^ 40±10 is allowed. 
The former implies that the light Higgs boson is ^ 110 GeV and, in principle, visi- 
ble to LEPII. The prediction for the Higgs boson mass in the tan /3 ~ 1 scenario is 
discussed and uncertainties related to (i) vacuum stability constraints, (m) differ- 
ent methods for calculating the Higgs boson mass, {Hi) two-loop calculations and 
{iv) GUT corrections are briefiy reviewed. It is shown that large left-right mixing 
between the i-scalars can significantly enhance the Higgs boson mass. That and 
an ambiguity in the size of the two-loop correction lead to our conservative upper 
bound of 110 GeV. Vacuum stability considerations constrain the i-scalar mixing 
and slightly diminish the upper bound (depending on the value of rrit). Improved 
two-loop calculations are also expected to strengthen the bound. 

Realizing the minimal supersymmetric extention of the standard model (MSSM) 
within a grand-unified theory (GUT) and assuming minimal matter content, often imply 
specific relations between the different Yukawa couplings. When considering the simplest 
example, i.e., "bottom-tau unification"^ hi, = h^, we find^ (in agreement with others"^) 
that the predicted ranges for as{Mz) and for tan /3 = Vupl^down are strongly correlated 

with the weak angle s'^[Mz) ~ 0.2324* and with mj, and are strongly constrained^. Such 
relations provide a consistency test of simple GUT models, as well as a means to probe the 
GUT structure. Furthermore, they can be used to distinguish such models from generic 
string models that typically do not constrain the Yukawa couplings in this manner. Thus, 
exploration of all implications of the bottom-tau unification (and similar relations) is well 
motivated. Below, we study the implications for the mass of the light MSSM (CP even) 
Higgs boson m/^o . We will motivate the exploration of the tan /3 — )• 1 limit and show that 
nifiO in this region is heavier than naive expectations, but most probably not too heavy to 
be seen in LEPII.'* Exploring that region, where the SM minimum oc cos^ 2/3 is shallow, 
requires us to pay special attention to vacuum stability constraints and we will briefiy 
comment on that issue. 

*Talk presented at SUSY-94, International Workshop on Supersymmetry and Unification of Fundamental 
Interactions, Ann Arbor, Michigan, May 14 -17, 1994. Pennsylvania Report No. UPR-0595T. 
*For ^ 143 GeV. 

ttts, tan /3 and mj are the strong coupling, the ratio of the two Higgs doublet expectation values and 
t-quark mass, respectively. 
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Fig. 1. The allowed region for tan /3 ~ 1. The i-quark pole mass m^"'^ is in GeV. The lower bound 
on tan (3 is from /i^ < 3 at the GUT scale, and the upper bound from bottom-tau unification 
(including theoretical uncertainties). The one s.d. range for m^°^^ ^ suggested by precision data, 
and tan /3 = 1 are indicated for comparison. Recent CDF analysis suggests m^"'^ ~ 174 ± 16 
GeV. For tan /3 < 1 no model with universal boundary conditions is possible. 

The most impressive evidence for TeV scale supersymmetry is the well publicized 
coupling constant unification.^ However, it is important to note that [up to matching 
functions induced by GUT-scale physics, which dominate the 0.008 theoretical error bar 
in (1)] unification requires as{Mz) ^ 0.12, i.e., 

as{Mz) = 0. 125 ±0.001 ± 0.008 + +3.2 x lO-^GeV"' [(mr'")' - (143 GeV)']. (1) 

The heavier the t-quark, the more the extracted value for the weak-angle is diminished and 
the predicted value for is increased. Two-loop contributions from Yukawa interactions 
Ha^ can diminish by as much as 0.003 (~ 2.5%) if the Yukawa coupling is of order unity 
(i.e., the function ^ —0.003). The above observations affect the region in parameter 
space which is consistent with the observational range for the 6-quark mass and with 
bottom-tau unification. In particular, for m^"'^ ^ 200 GeV only two regions are allowed: 
(i) tan /3 1 (see Fig. 1) and (ii) large tan /3. For example, if mf°'^ 174 GeV^ then 
tan /3 1.5 or tan /3 sa 40 ± 10. Hereafter, I will discuss region (z) only. An important 
implication of Yukawa unification, choosing the tan /3 near unity region, is that the SM- 
like Higgs boson mass nearly vanishes at tree-level. Thus, a careful examination of the 
loop corrections is required. 

There are important theortical uncertainties which constrain the predictive power 
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Fig. 2. The correction to the 5-quark MS mass mfj{Mz) ~ 3 GeV (in GeV) from finite 
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174 GeV and universal boundary conditions for the soft param- 



eters at the unification scale are choosen at random. All points included are consistent with 
electroweak breaking, tan /3 = 20 is indicated for comparison. For tan /3 < 20 corrections are less 
than 10%. 



of the model. ^ These are taken into account (assuming no conspiracies), e.g., in (1) and 
Fig. 1, but will not be discussed here. The following points, however, should be stressed: 

• cts dependence: We use the predicted value of as(Af^)[s^, mt, ...], e.g., as(Af^) 0.12 

for mf' ^ 100 GeV and a,{Mz) ^ 0.13 for mf'" ^ 180 GeV. Using instead a fixed 
value a°(Af^) for the whole nit range implies some hidden assumptions on GUT- 
scale corrections. The lower the value of (and the heavier the t-quark) the larger 
and less likely are these corrections. 

• Finite-loop corrections to mb[MsusY)- The relevance of those corrections was pointed 
out recently. However, while being extremely important for large tan /3 (maybe even 
to the point of removing any predictive power), they are negligible in region (z), as 
can be seen in Fig. 2. The details of the low-energy spectrum are nearly irrelevant 
in determining the size of that region. 

Different aspects of region (z) were discussed by various authors. The possibility 
that the Higgs boson mass is induced only at the loop-level was discussed by Diaz and 
Haber.® The consistency of that scenario with bottom-tau unification recently led to 
greater attention to that region of parameter space. "^'^'^^ The consistency with Yukawa 
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unification is easily understood: From perturbativity one has a mj-dependent lower bound 
on tan /3 ^ 1. Thus, we are in a region of large top- Yukawa coupling hf ~ f. Slightly 
decreasing tan /3 (i.e., decreasing sin /3 or increasing hf) at Mz would lead to divergences 
(i.e., hf "hits" its Landau pole) at higher scales^^: hf has a quasi-fixed point at tan /3 f . 
This balances too large corrections to the hb/h-j- ratio, can explain the heavy t-quark, 
and makes region (z) relatively insensitive to theoretical uncertainties. 

For tan /3 = f there is a custodial symmetry SU{2)l x SU{2)r — )• SU{2)l^r in 
the Higgs sector. The deviation of tan /3 from unity measures the symmetry breaking, 
which is induced at the loop level due to the Yukawa interactions. Thus, the enhanced 
symmetry can be responsible for the smallness of the region. It also puts an upper bound 
on mf°'^ ^ 185 GeV so that tan (3 < 2. 

As mentioned above, we will not discuss all possible signatures, but will focus on the 
loop-induced Higgs boson mass. (Another interesting signature is a possible light mixed t- 
scalar.^°) In particular, Barger et al. argued, ^"^ using the leading-logarithm approximation, 
that m/jO ^ 85 GeV for ° ^ ^ 160 GeV. We will show that this is rather a typical mass 
range and that the upper bound, though still relevant for LFPII, is much higher. An 
important enhancement comes from the large left-right mixing in the t-scalar sector and 
from the large split between the two t-scalar mass eigenstates. That enhancement leads 
to the upper bound m/^o ^ 100 (110) GeV for m^"'^ ^ 160 (175) GeV. The enhancement 
is sensitive to the type of vacuum stability constraints one imposes (see below). 

From the minimization of the Higgs scalar potential one has for the supersymmetric 
Higgsino mass parameter //^ oc l/[tan^/3 — 1], and //^ diverges as tan /3 — )• 1, i.e., in the 
symmetric limit. (In fact, one expects finite loop corrections to be relevant near the limit, 
and we can assume, without loss of generality, tan /3 ^ 1.1.^"*) The large //-parameter 
dictates the phenomenology of the scenario. In particular, the tree-level CP even mass 
matrix is nearly degenerate: it has a very heavy and a nearly zero mass eigenvalue (this is a 
trivial consequence of the custodial symmetry). In practice, the light tree-level eigenvalue 

with the lower bound on tan/3). 
. The upper bound on the loop 



m 



< Af^|cos2/3| is small but grows with m^"'^ (i.e 



At the one-loop level, one has m\o = mjo'^ + 



A2 



correction can be estimated^^ (assuming only t-scalar contributions) 



m 



< M|cos2 2/3 + 



47rs2 1 



0M| 



In 



where 



A. 



-|- m 



m 



t2 



mi 



m 



t2 



4m? 



^iff^ln 



In 



(3) 



Due to the large // parameter the left-right mixing 6t between the two t-scalars can be 
substantial and the leading-logarithm and the mixing Ag^ terms in (2) can be equivalent. 
Thus, m/jO is enhanced by an additional factor of ~ \/2. However, since tan j3 increases 
with nit (from perturbativity), and Ag^ decrease and we obtain an interesting interplay 
between the overall factor of and Ag^ . 

Fq. (2) is, of course, only an approximation. The results of a complete calculation 
for the Higgs boson mass are shown in Fig. 3 (for the details of the calculation, see Ref. 
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Fig. 3. The Higgs boson mass nif^o (in GeV) distribution in a sample of monte-carlo calculations 
with random universal boundary conditions at the unification scale and using the effective po- 
tential method with subtraction scale Mz- tan/3 is constrained as in Fig. f. m^"'^ is in the range 
(a) [f55, f65], (5) [f65, f75], (c) [f75, f85] GeV and has the gaussian distribution f43± f8 GeV. 



4). tan /3 is constrained as in Fig. 1, we assume the gaussian distribution ^ = 143 it 18 

GeV (from precision data+) and all superpartner masses are constrained to be ^ 1 TeV. 
The monte-carlo generated histograms contain information on the upper bounds and on 
the distributions. The upper bound is a function of nit'. 

• niho ^ 100 GeV for mf'" ^ 160 GeV, 

• niho ^ 110 GeV for mf' ^ 175 GeV, 

• niho ^ 120 GeV for mf'" ^ 185 GeV. 

The distribution has two peaks, one which is enhanced by the t-scalar mixing term and 
which determines the upper bound, and a peak at a much lighter mass from points with 

+ Using the recent CDF range'^ for would only change the relative population in the different 

histograms and would not affect our discussion. 
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Fig. 4. The upper bound on the Higgs boson (MS) mass m/^o as a function of ^ (all masses 
are in GeV). In the constrained case tan [i is constrained by bottom-tau unification as in Fig. 1. 
The general case is for all tan/3 (consistent with electroweak breaking). Note the local minimum 
in the constrained curve which is due to the interplay between the mixing enhancement and 
the m'l factor. The upper bound curves are derived using monte-carlo methods and are not 
rigorous. Generic upper bounds in the MSSM and in non-minimal models (NMSSM) derived 
from perturbativity considerations are shown for comparison. Also shown a suggestive lower 
bound in the non supersymmetric case (we use = 0.125 in the formula of Ref. 17). 

no enhancement. Thus, the mass is typically much lighter than the upper bound, and, if 
the t-quark is lighter than ~ 170 GeV, may still be relevant for LEPI. (Note that we do 
not attempt to give rigorous bounds.) 

We compare the above bounds to those derived for any tan fi in Fig. 4. The ~ 30 GeV 
difference is due to the non- vanishing tree-level mass ^ Mz in the general case. (Loop 
corrections are, however, typically smaller in the general case.) One can also compare 
either case with typical upper bounds m/^o ^ 130 — 150 GeV derived from perturbativity 
considerations,^^ and to the lower bound given by Sher^ for the non-supersymmetric case, 
m/jO ^ 132 GeV. However, strong assumptions regarding vacuum stability are made in 
the latter case. 

Complications in the calculations, which are explored elsewhere,'* suggest even stronger 
upper bounds: 

1. The magnitude of the two-loop terms is known to be small and probably neg- 
ative.*® However, the exact magnitude is ambiguous given the complicated low-energy 
structure of the model. For example. In Fig. 3 we used the effective potential method*^ 
to extract m/^o. In that method A/^o is given to a one-loop leading-logarithm order. It is 
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straightforward to show that the renormahzation-point dependence of the one-loop lead- 
ing logarithm [e.g., in (2)] is equivalent to two-loop next-to-leading logarithms which are 
typically positive. Thus, one could overestimate m° when choosing the renormalization 
point Q = Mz (as in Fig. 3). In Fig. 5a Q = 600 GeV, and indeed lower masses are 
suggested. In particular, if the scalar quarks are of the order of TeV or heavier, one would 
falsely get large Higgs boson masses using that method, unless, Q is properly adjusted 
(so that the loop expansion does not break down). Alternatively, we could use different 
methods, e.g., the renormalization group method^° leads to typically lower one-loop upper 
bounds (in agreement with Fig. 5a). This is because one is implicitly including negative 
two-loop leading logarithm contributions in that method. In short. Fig. 3 corresponds 
to a conservative estimate of the upper bounds. Once the non-trivial task of a complete 
two-loop calculation is carried out, we expect slightly stronger upper bounds. 

2. In some cases large t-scalar mixing corresponds to the SM model populating only 
a local minimum of the full scalar potential. That is, when considering all scalar fields 
on equal footing, the global minimum may be in a direction in which the t-scalar has a 
vacuum expectation value and, thus, the physical vacuum is not stable. Fliminating cases 
in which a global color and charge breaking minimum exists slightly reduces the upper 
bound and also affects the mixing peak population. The effect is reduced with increasing 
nit and never diminishes the mixing peak completely. One could try to identify dangerous 
directions in the multidimensional scalar-field space and characterize those directions by 
analytic constraints. However, such constraints entail assumptions that are not always 
useful for identifying the deepest minimum (i.e., the constraints are too weak). In addition, 
analytic constraints typically do not distinguish a local minimum from a global one and 
can exclude legitimate points (i.e., are too strong). Thus, analytic treatments are often not 
satisfactory. Furthermore, a point in parameter space which is consistent with electroweak 
symmetry breaking contains at least a local minimum which does not respect color and 
charge symmetries. This is shown in Ref. 4 where we also compare numerical and analytic 
treatments of the problem. Recent progress was also reported by Carena and Wagner. 

Lastly, the prediction of m/^o is relatively stable when considering model-dependent 
evolution of the soft parameters between the Planck and GUT scales. This is shown in 
Fig. 5 and is due to the cancellations between corrections to the relevant soft parameters 
and to /i^. (However, the t-scalar mixing and mass are affected.) Thus, our discussion is 
only slightly affected by that uncertainty. 

To conclude, we studied the one-loop Higgs boson mass in the region of parameter 
space which is well motivated by Yukawa unification. We treat all theoretical uncertainties 
in a reasonable manner and the upper bounds given, e.g., m^o ^ no GeV (for ^ 175 

GeV), correspond to the most conservative estimate. The Higgs boson in this scenario is 
most likely in a range that is, in principle, visible to LFPII [unless the t-quark is heavier 
than ~ 180 GeV, in which case there is not much motivation to consider region (z)]. Two- 
loop calculations and vacuum stability constraints slightly strengthen the upper bounds. 
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Fig. 5. The Higgs boson {mf^o) vs. Higgsino {fi) mass plane (in GeV) for = 174 GeV, 

tan/3 = 1.52, and with random universal boundary conditions at (a) the GUT scale ~ 10^^ 
GeV, (b) at the (reduced) Planck scale ~ 10^® GeV. In case (b) we assume a minimal SU5 model 
above the GUT scale (for details, see Ref. 22). The effective potential method with a subtraction 
scale Q = 600 GeV is employed, which leads to lower values of nif^o than for Q = Mz- 
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